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ABSTRACT 

We consider the anchor placement problem in localization 
based on one-way ranging, in which either the sensor or the 
anchors send the ranging signals. The number of anchors de- 
ployed over a geographical area is generally sparse, and we 
show that the anchor placement can be formulated as the de- 
sign of a sparse selection vector. Interestingly, the case in 
which the anchors send the ranging signals, results in a joint 
ranging energy optimization and anchor placement problem. 
We make abstraction of the localization algorithm and instead 
use the Cramer-Rao lower bound (CRB) as the performance 
constraint. The anchor placement problem is formulated as 
an elegant convex optimization problem which can be solved 
efficiently. 

Index Terms — Anchor placement, ranging energy opti- 
mization, localization, sparsity, convex optimization. 

1. INTRODUCTION 

Localization is an important and extensively studied topic in 
wireless sensor networks (WSNs) [1]. Localization can be 
performed using a plethora of algorithms [1] (and references 
therein), which exploit inter-node measurements like time-of- 
arrival (TOA), time-difference-of-arrival (TDOA), angle-of- 
arrival (AOA), or received signal strength (RSS). 

The performance of any location estimator depends not 
only on the algorithm but also on the placement of the an- 
chors (nodes with known locations). Anchor placement is a 
key challenge in localization system design, as certain anchor 
positions not only deteriorate the performance but also result 
in ambiguities or identifiability issues. In [2], the effect of 
anchor placement is studied using the geometric dilution of 
precision (GDOP). The idea of GDOP is borrowed from the 
global positioning system (GPS) and it is obtained from the 
Cramer-Rao lower bound (CRB) with simplifying assump- 
tions on the noise model, i.e., equal variances on all the range 
estimates. This assumption is valid for GPS due to the ap- 
proximately equal distances between the anchors (the satel- 



lites) and the nodes, but it does not translate well to WSNs 
where the noise variance is proportional to the distance (typ- 
ically different) between the anchors and the sensor. In [3], 
the effect of the anchor positions has been studied by empiri- 
cally identifying the ambiguity issues. However, they do not 
provide any algorithm for anchor placement. 

The anchor placement problem can be interpreted as the 
problem where we divide a specific anchor area A into M grid 
points and select the positions of the K anchors as the best K 
grid points out of M grid points, where K <C M. Here, the 
K selected anchors are deemed the best, if they guarantee a 
certain minimal accuracy on the location estimates within a 
specific sensor area S. In practice, K is not known, and this 
makes anchor selection a combinatorial problem involving an 
exhaustive search over all the 2 M possible anchor positions, 
and the computation over S is even more cumbersome. 

In this paper, we consider TOA-based ranging, however, 
we do not restrict ourselves to a particular localization algo- 
rithm. Instead, we use the CRB as a performance constraint. 
The anchor placement problem is studied for the following 
two cases of TOA-based one-way ranging: a) the anchors 
send the ranging signals (OW-A) and b) the sensor sends the 
ranging signals (OW-S). The anchor placement problem is 
formulated as the design of a sparse selection vector. For 
the OW-A case, the sparse solution yields the ranging ener- 
gies 1 that the anchors should adopt leading to a solution for 
the joint ranging energy optimization and anchor placement 
problem, and for the OW-S case the sparse solution yields the 
optimal anchor positions. Using the CRB as a performance 
constraint, we can formulate the anchor placement problem 
as a semidefinite programming (SDP) problem. The major 
contribution of this paper is the novel framework for anchor 
placement in WSNs exploiting the inherent sparse nature of 
the spatially distributed anchors. 

2. SYSTEM MODEL AND PRELIMINARIES 

We consider a two-dimensional network with one sensor lo- 
cated in the sensor area S and M possible anchors located at 
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'More details on ranging energy optimization can be found in [4], how- 
ever, the design of optimal anchor positions is not considered in [4]. 



the M grid points of the anchor area A. Let the coordinates 
of the sensor and the mth anchor be denoted by the 2x1 
vectors s and a m , respectively, where s is assumed to be un- 
known but known to be within S. We further assume that all 
the nodes are time-synchronized (using [5] for instance). Let 
the pairwise distance between the sensor and the mth anchor 
be denoted by d(a m , s) = ||a m — s|| 2 . In practice, the pair- 
wise distances are obtained by ranging and they are generally 
noisy. For the sake of simplicity, we consider a TOA-based 
one-way ranging. 

2.1. Anchors send the ranging signals: OW-A 

The ranging signals sent by the anchors are scheduled such 
that they can be separated at the sensors. Let the mth anchor 
broadcast a ranging signal ^fe^s a {t) of energy e m at time 
T m>8 , and upon reception at the sensor, the TOA T s>m is esti- 
mated. The range estimate of the sensor to the mth anchor is 
then given by 



i(3m|S) — c(T m , s T s 



(1) 



where c is the propagation speed of a wave in the medium. 
Under the assumption that a line of sight (LOS) channel ex- 
ists between the sensor and the mth anchor, it is motivated 
in [4,6] that the range estimate d a (a m ,s) is Gaussian dis- 
tributed as A/"(d(a m , s), a 2 n s ) with variance cr^ 



-l£a 
e m 7 2 ■ 



Here. 



Pa = — pT- With Fl = 



/~ oo (27r_F) 2 |Sa(-F)| 2 d-F 



the 



mean square bandwidth of the ranging signal (S a (F) being 
the Fourier transform of s a (t)), and with N s /2 the two-sided 
power spectral density (PSD) of the additive white Gaussian 
noise (AWGN) at the sensor. Further, the signal suffers an at- 
tenuation of 7 2 = ad(a m , s)~@, with a and f3 the path gain 
at 1 m and path-loss coefficient, respectively. 

2.2. Sensor sends the ranging signal: OW-S 

In the OW-S case, the sensor broadcasts a ranging signal 
Y / eJs(t) of energy e s at time T s>m , and upon reception at the 
mth anchor, the TOA T m ,s is estimated. As earlier, the range 
estimate of the sensor to the mth anchor is given by 



d s (a m ,s) = c(T s , 



T„ 



(2) 



which is again assumed to be Gaussian distributed with 



variance at 



e s where p s = c N - a J 2 with F 2 



IZ,(^ F Y\s(F)\ 2 dF 

— joo \s(F)\ 2 dF — mean square bandwidth of the rang- 
ing signal (S(F) being the Fourier transform of s(t)), and 
with N a /2 the two-sided PSD of the AWGN at the anchor. 

3. PERFORMANCE MEASURE 

We make abstraction of the localization algorithm, however, 
we assume that the TOA estimates are unbiased and achieve 
the CRB asymptotically. Even though this assumption is too 



optimistic for a practical system, the CRB has a very attrac- 
tive mathematical structure resulting in a selection problem 
that can be efficiently solved using convex optimization tech- 
niques. Moreover, the CRB optimal anchor positions can im- 
prove the performance of any practical localization algorithm. 

The CRB for OW-A and OW-S can be derived based on 
the CRB for TOA-based two-way ranging [4, 6]. 

3.1. CRB for OW-A 

For OW-A, the CRB of s is denoted by C a £ R 2x2 , and is 
computed as follows 

M 

C- 1 = F a (e, s)=J2 e m F a , m (s) (3) 

m— 1 

where F a , m (s) = ap~ 1 d{a m , s)~^~ 2 (s - a TO )(s - a m ) T 
and F a (e, s) is the Fisher information matrix (FIM). Here, the 
vector e = [ei, eg, . . . , eNi\ T is the joint selection and anchor 
ranging energy vector that has to be designed, where a non- 
zero entry of e not only indicates that the anchor position is 
selected but also represents the ranging signal energy that the 
selected anchor should adopt. 

3.2. CRB for OW-S 

Similarly for OW-S, the CRB of s is denoted by C s el 2x2 , 
and is computed as follows 



M 



C; 1 = F s (w,s) = J2 w m F s , m {s) 



(4) 



where F s>m (s) = e s ap a 1 d(a m ,s) 13 2 (s - a m )(s - a m ) T 
and F s (w, s) is the FIM. Here, w = [wi, . . . ,wm] T € 
{0, 1} M is the selection vector to be designed, where w m = 
1 (0) indicates that the related anchor is (not) selected. 
3.3. Identiflability and ambiguity 

The FIM is singular when both the anchors and the sensor 
are collinear for which the CRB will be infinity. Hence, a 
CRB -optimal anchor placement would avoid the CRB being 
infinity. However, when only the anchors (even if M > 3) 
are collinear (excluding the sensor) a "mirror" ambiguity is 
obtained, and this effect cannot be seen with the FIM, and sur- 
prisingly the FIM will be non-singular. Hence, a (local) CRB- 
optimal anchor placement can still result in such ambiguities. 
Such solutions can be avoided with additional constraints or 
prior information on the parameters, e.g., by a constraint on 
the sensor area, such as an orthant or a half plane. 

A constrained CRB generally gives a lower bound for pa- 
rameters with such deterministic constraints. A constrained 
CRB is derived from the unconstrained CRB and a non- 
redundant constraint set. However, for orthant or half plane 
constraints, it is shown in [7] that the unconstrained CRB 
is the same as the constrained CRB, and does not yield a 
lower bound. Hence, throughout this paper, we use the un- 
constrained CRB in (3)-(4) as a performance measure and we 
will select our S carefully. 



3.4. Constraint for accurate positioning 

Using the definition of accurate positioning from the Federal 
communication commission (FCC) [1], for every s within S 
we constrain the localization error £ = s — s to be within an 
origin-centered circle of radius R e with a probability higher 
than P e , i.e., Vs e S,Pr(||£|| 2 < R e ) > P e . The values of 
R e and P e define the accuracy required and are assumed to be 
known. 

Sufficient conditions to satisfy this accuracy requirement 
for OW-A and OW-S are X min (e, s) > A and A OT j„(w, s) > 
A, respectively [1,4]. Here, A m j n (e,s) and A TO j„(w, s) are 
the smallest eigenvalues of the matrices F a and F s , respec- 
tively, and the threshold A for a Gaussian distributed location 
estimate is given by A = ln( ^ p ) and for an unknown 

distribution by A = JjCjzrp-) [1.4]. 

More specifically, we lower bound the eigenvalues of 
the FIMs F (e, s) and F s (w, s), so that the related CRB is 
upper bounded. The inequality constraints A m j„(e, s) > A 
and A m i„(w,s) > A are equivalent to the following lin- 
ear matrix inequalities (LMIs): Yl m =i e mFa,m( s ) t AI2 
and J2 m =l w m~Fs,m{s) h AI2, respectively, where A >; B 
means that A — B is a positive semidefinite matrix, and 1 2 is 
a 2 x 2 identity matrix. It is well known that the solution set 
of e and w satisfying these respective LMIs is convex [8]. 

4. SPARSITY-EXPLOITING ANCHOR PLACEMENT 

Let us now assume that M possible anchors are placed on a 
discrete grid obtained by uniformly sampling an anchor area 
A. And remember that the location of the sensor is unknown 
within a sensor area S. In many localization applications, it 
makes sense to assume that a limited number of anchors ser- 
vice a prescribed geographical area. This assumption natu- 
rally leads to the design of a sparse vector for optimal anchor 
placement. More specifically, we aim to design a sparse joint 
selection and ranging energy vector e and a sparse selection 
vector w for OW-A and OW-S, respectively. 

4.1. OW-A: Joint anchor placement and ranging energy 
optimization 

When the anchors send the ranging signals, the optimization 
problem can be written as 



mm e 

66R M 

M 



(] 



s.t. e m F a , m (s) - AI 2 >r 0, Vs e S 

m— 1 

OmxI < e < CblA/xi, 



(5) 



where the £o-(qaasi) norm refers to the number of non-zero 
entries of e, i.e., ||e|| := \{m : e rn 7^ 0}|. In addition to the 
performance constraint, an energy source is positive-valued, 
and is generally constrained by a prescribed value et due to 
the practical limitations of a source. Here, the threshold A 



is indirectly the sparsity-inducing parameter (for a fixed e^), 
where implies a sparser solution. 

It is well known that the l^-norni optimization is NP-hard 
and non-convex. A computationally tractable solution is to 
use the traditional convex surrogate for the ^-novm, which 
leads to the following optimization problem: 



e* = are min 1 e = 1 1 e I L 

eel 1 ' 1 

M 

s.t. 

a,m 

(s) - AI 2 h0, Vs g S 

m— 1 

Omxi < e < e&lMxi- 



(6) 



The purpose of the ^i-norm in (6) is twofold: promote the 
sparsity of the spatially distributed anchors and more impor- 
tantly minimize the "total" ranging energy of the network. 
This is a SDP problem in a standard dual form, and its solu- 
tion upper bounds the dual feasible e, i.e., e* < l T e. 

The optimization problem in (6) has multiple global so- 
lutions, due to the fact that the £i-norm is not strictly con- 
vex. Even though the £i-norm minimizes the total ranging 
energy of the network, it does not necessarily minimize the 
number of anchor nodes. On the other hand, the non-convex 
(intractable) £o- noim optimization leads to a higher energy. 
To improve upon the £i-norm solution, we propose an alter- 
native (convex) optimization algorithm, which also results in 
a correct solution, but with a fewer number of anchors. For 
this purpose, we modify (6) and use the sparsity-enhancing it- 
erative re-weighted ^i-norm minimization [9] originally used 
for linear t\ -regularized least-squares (LS) problems in com- 
pressed sensing (CS). 

4.2. Sparsity-enhancing iterative algorithm for OW-A 



The iterative re-weighted £i-norm minimization [9] algorithm 
is adapted to suit our problem. Let u = [ui, . . . , um] t €E K m 
denote the weight vector. The iterative algorithm proceeds as 
follows: 

1. Initialize the iteration counter k — and the weight vector 
u<°> = 1m. 

2. Solve the weighted ^i-norm minimization problem 



(k)T (k) 

mm xv ' e y ' 

M 

s.t. e^ ) Fa, m (s) - AI 2 >r 0, Vs e 5 



Oa 



< e v < e 6 ljv/xi 



(7) 



for each i 



for the optimum e' fc ' in the fc-th iteration. 

(k) 1 
Update the weight vector u\ — r 

e +l e i 1 

1,...,M. 

Stop on convergence, or when k attains a specified maximum 
number of iterations k max , else, increment k and go to step 2. 
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Fig. 1: Joint anchor placement and energy optimization for OW-A, for P e — 0.95, et, = 10J, and R e = 4cm. 



The weight updates force the small entries of the vector 
to zero and avoid inappropriate suppression of larger entries. 
The parameter e > is the threshold which provides stability, 
and guarantees that the zero valued entry of |e( fe )| does not 
strictly prohibit a nonzero estimate at the next step. 

4.3. OW-S: Anchor grid points selection 

In the OW-S case, where the sensor sends the ranging signals, 
a similar optimization problem can be formulated for the se- 
lection of the anchor grid points, which is given by 



l 1 w = llwl 



f b = arg mm 

M 

s.t. w m F s . m (s) - AI a hO, Vse S. 



rn—1 



(8) 



The optimization problem in (8) is a non-convex Boolean 
programming problem. However, this can be brought to 
the standard dual SDP form using the Lagrangian relax- 
ation w m (w m — 1) = 0, and introducing a new variable 
W = ww T with elements [W] mn ,m,n = 1,...,M [8]. 
The optimization problem is then given by 



(w* dp , W* d ) =arg min 1 w=|l w lli 

wa M ! w6i MXM 



s.t. ™mF s , m (s) - AI 2 h 0, Vs G S 

m=l 

[W] mm = w mj m = l,...,M, 
W wl^ _ 



(9) 



Here, the rank-1 constraint on W is relaxed, and (9) can be 
solved efficiently in polynomial time. The SDP in (9) pro- 
vides a good approximation for the Boolean problem in (8), 
and the solutions for (8) and (9) are upper bounds for the dual 



feasible w, i.e., w^ p < 



sdp 



< 1 w. From w 



sdp 



the ap- 



proximate Boolean solution to w can be obtained using ran- 
domization techniques. 



4.4. Sparsity-enhancing iterative algorithm for OW-S 

The sparsity-enhancing iterative algorithm for OW-S can be 
derived following similar lines as discussed in Section 4.2, 
and it proceeds as follows: 

1. Initialize the iteration counter k = and the weight vector 
u<°> = 1m. 

2. Solve the weighted i!i-norm minimization problem 



(k)T 

mm u v ' w 

wei M ,W6« M ^' 




M 




S.t. ^2 WmFa.rr 


(s) - AI 2 >r 0, Vs e S 






[W] mm = w 


Tl , m = l,...,M, 




W w " 

w T 1 


y o. 



(10) 



for the optimum w« in the fc-th iteration. 



(*■) 



for each i = 



3. Update the weight vector u o 
1,...,M. 

4. Stop on convergence, or when k attains a specified maximum 
number of iterations k max , else, increment k and go to step 2. 

The Boolean solution can again be obtained using randomiza- 
tion. 

Note that the sensor ranging energy e s is not optimized. 
Once the optimal anchor grid points are obtained, then e s can 
be easily optimized on the reduced size problem. 

5. SIMULATION RESULTS 

To test the proposed SDP-based algorithms, we use CVX [10]. 
CVX in turn calls SeDuMi, a MATLAB implementation of 
the second-order interior-point methods for computations . 

We simulate a network with M — 80 anchor grid points 
and for two scenarios of the sensor area S consisting of 608 
and 25 sensor grid points as shown in Fig. 1 and Fig. 2, re- 
spectively. Here, we have sampled S to generate discrete sen- 
sor grid points 2 . We use the following parameters for the sim- 
ulations: a = 1, = 2, c = 3x 10 8 ms _1 , = 8 GHz, 



2 A more profound treatment of the gridding of the sensor area can be 
found in [4]. 
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domized). pig 2: Anchor placement for OW-S, for P e = 0.95, e s = 10J, and R e = 5cm. 



N s /2 = OdBW/Hz, and e = 1(T 8 . The constraint on the 
anchor ranging energy is e& = 10J, and the sensor ranging 
energy e s = 10J is used. 

The optimal anchor placement for the OW-A case, based 
on ^-norm minimization of (6) and the iterative ^-norm 
minimization is illustrated in Fig. lb and Fig. lc, respec- 
tively. The algorithm selects the anchor grid points close to 
the sensor area due to the assumed path-loss model. There 
exits many similar solutions, however, we are interested in 
the best solution. This effect can be seen in the not so sparse 
£i-norm solution, while the iterative £i-norm solution en- 
hances the sparsity, yet yielding a correct solution. Fig. la 
shows the optimal ranging energies that the anchors should 
adopt. 

Fig. 2a illustrates the anchor selection for a circular an- 
chor grid with M = 80 points for OW-S. Here, we con- 
sider a relatively small sensor area. For the £i-norm solu- 
tion, all the grid points are nominally the same, and the se- 
lection weights are spread over all the nodes. However, the 
iterative £i-norm solution enhances the sparsity. Fig. 2a also 
shows an approximate Boolean solution computed using the 
solution from the iterative algorithm followed by randomiza- 
tion. Note that the randomization using the ^i-norm solution 
will not yield a meaningful Boolean solution, as the £i-norm 
solution is not sparse for this scenario. Fig. 2b shows the 
corresponding anchor placement using the iterative ^i-norm 
algorithm followed by randomization. 

An exhaustive search based algorithm to select 14 grid 
points out of 80 for the case shown in Fig. lc would need 
( 8 g) x 608 w 10 17 searches over the constraint set, which 
is clearly intractable. On the other hand, the complexity of 
the proposed algorithms for both OW-A and OW-S cases are 
polynomial in the number of variables and the number of con- 
straints. 
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